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Abstract 

We determine the statistical properties of wave functions in disordered quan- 
tum systems by exact diagonalization of one-, two- and quasi-one dimensional 
tight-binding Hamiltonians. In the quasi-one dimensional case we find that 
the tails of the distribution of wave-function amplitudes are described by the 
non-linear er-model. In two dimensions, the tails of the distribution function 
are consistent with a recent prediction based on a direct optimal fluctuation 
method. 
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It is well established that disordered quantum systems in the metallic regime (i.e., in 
the limit of weak disorder) and highly excited classically chaotic quantum systems exhibit 
universal quantum fluctuations that can be described by random matrix theory (RMT): 
statistical properties, on the scale of the mean level spacing, of eigenvalues, eigenf unctions, 
and matrix elements are universal, i.e., they do not depend on the microscopic details of the 
systems under consideration [l|-[5|]. 

However, in ballistic, classically chaotic quantum systems, non-hyperbolic phase-space 
structures may lead to deviations from universal RMT statistics ||. Similarly fluctuations in 
disordered, classically diffusive quantum systems may deviate considerably from the RMT 
predictions due to increased localization. This effect is naturally very significant in the 
tails of distribution functions |J (corresponding to rare events) of wave-function amplitudes 
|7HT3|, of the local density of states [|7|,|T2|] , of inverse participation ratios fi"2| , |T3|l and of 



NMR line shapes 0. In all of these cases (with the exception of Ref. Jj]] which deals with 
one-dimensional (ID) systems), the distribution functions have been calculated using the 
non- linear cx-model (NLSM). Very recently, this approach has been extended to ballistic 
systems [0,0] (see also [|0-[T§). 



In Ref. |TP[] a direct optimal fluctuation method |2D| was used to calculate the tails 



of distributions of current relaxation times and wave-function amplitudes; and predictions 
differing from |p Hl^1 were put forward. This led the authors of [19|] to question the suitability 
of the NLSM to describe rare events in disordered conductors. 

It is thus of great interest to test the predictions of |7|-[T3| and ]19|] against results of 



independent calculations. In this letter, we have determined distribution functions of wave- 
function amplitudes by exact diagonalization of ID, 2D and quasi-lD tight-binding Hamil- 
tonians; in this case rare events correspond to unusually high splashes of wave-function 
amplitudes. We note that wave-function amplitude distributions can be measured in micro- 



wave experiments 21,22 



We use the Anderson model of localization K3] which is a tight-binding model on a 
(i-dimensional hyper-cubic lattice 



H = t rr' C l C r > + Yl V rC ] r C r . (1) 

Here r — (x,y, . . .) denotes sites on the lattice, cj. and c r are the usual creation and anni- 
hilation operators, the hopping amplitudes are t rr r = 1 for nearest neighbour sites and zero 
otherwise. The on-site potential v r is taken to be uncorrelated white noise, with zero mean 
and variance (tyiv) — 5 rr 'W 2 /12. The parameter W characterizes the disorder strength. 
As is well-known (see for instance [2^25|), the eigenvalues Ej and eigenfunctions ipj(r) of 
this Hamiltonian, in the metallic regime, exhibit fluctuations described by RMT. In this case, 
Dyson's Gaussian orthogonal ensemble Q is appropriate. When the matrix elements t rr > 
of (|l]) are given an appropriate complex phase factor, Dyson's unitary ensemble Q applies. 
We refer to these two cases by assigning, as usual, the parameter f3 = 1 to the former and 
(3 — 2 to the latter. The metallic regime is characterized by g ^> 1 where g = IkvV DL~ 2 
is the dimensionless conductance (we take fi = 1). Here v = 1/(VA), A is the mean level 
spacing and V the volume. D = v\rjd is the diffusion constant, r the mean free time and 
vp the Fermi velocity. Four length scales are important: the lattice spacing a, the linear 
extension L, the localization length £ and the mean free path i = v-pr. 

By diagonalizing the Hamiltonian H using a modified Lanczos algorithm ||26|| , we have 
determined the distribution function 

ME,r ] t) = A(Y J S(t-\^(r)\ 2 V)g v (E-E,)) . (2) 
\ — / w 
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Here (■ ■ -)w denotes an average over disorder realisations. The wave functions are normalized 
so that (|V'i( r ')| 2 )vF = V~ l an d g^E) is a window function of width rj, centered around E — 
and normalized to unity. In the following we describe the results of our calculations and 
compare them to the predictions of Refs. [pHT^,P^ • 



ID case. The eigenstates in a disordered chain are localized with localization length 
£ = 4£. According to Ref. |7j the distribution of wave-function amplitudes in a disordered 
chain of length L is [^7| 

/(*;*) =4 «p(-$ (3) 



for L ^> £, independent of x and /3. Our results for (f(E,x;t)) x ((• • -) x denotes an average 
over the x coordinate) in Fig. [1] show very good agreement with Eq. (^|) for large L/C,. 
The deviations at small t for L/£ = 4.76 are due to the fact that Eq. (|3]) is only valid 
asymptotically for large L. It does not take into account that in a system of finite length 
L, the smallest amplitude of a normalized, exponentially decaying wave function is of the 
order of t c ~ {L/C) exp(— L/C,). This cut-off is shown in Fig. [l|. 

Quasi- ID case. In this case, as was shown in Refs. f§||, the NLSM can be solved exactly 
for the distribution function fp{E, x; t), using a transfer matrix approach p8 |. The result is 

^^^^r^ y ^ + ^- (4) 

f 2 (E,x ] t) = ^- 2 Y(t). (5) 
Here Y (z) =W(z£/L,x/£)W(z£/L,(L — x)/£) and W(z,r) obeys the differential equation 

with initial condition W(z, 0) = 1. The function W(z, t) may be determined in terms of 
an eigenfunction expansion of the operator z 2 d 2 z — z P, |i~2| , |i~3| . The change of the body and 
the tails of distribution fp{E, x; t) due to increasing localization may thus be parameterized 
by a single parameter which we define to be X = (f3/2)L/£ j . Here £ = fiwvDS where S is 
the cross-section of the wire. Thus X does not depend on (3. In the metallic regime (where 
X -»• 0) Y(z) ~ exp(-z) which leads to the usual RMT results A (0) (t) = exp(-t/2)/v / 2vrt 
and (t) = exp(— t). The former distribution (/? = 1) is often referred to as the Porter- 
Thomas distribution M. For increasing localization (finite but still small X), the 0(X)- 
corrections to the body of the distribution function fg are obtained by expanding Y ~ 
exp(— t)[l + (3~ 1 t 2 P(x, x; 0)] where P(x,x';u;) is the one-dimensional diffusion propagator. 
The result is f p (E,x;t) = /^ 0) + 8fp{E,x;t)} with 
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8fp{E,x;t)~P{x,x;0)- 



3/4-3t/2+t 2 /4 tor 13 = 1, 
1 — 2t+t 2 /2 for p = 2, 

valid for t < X~ 1/2 . In the tails (i > X' 1 > 1) of f (E,x;t), Eqs. (§ and (D simplify to 
1 

/>(£, x; t) ~ A p (x, X) exp(-2f3^t/X) . (7) 



This result may also be obtained within a saddle-point approximation to the NLSM |L0 
The prefactors Ap(x, X) for j3 = 1, 2 are given in p|JT0[ . 



According to Refs. |||28] and the NLSM applies provided the following conditions 

are satisfied: 

K M < fep-S < A; F L (8) 

where fep is the Fermi wave vector and 5 is the cross section of the wire. The first condition 
ensures that disorder is sufficiently weak. The second condition implies that apart from the 



sample geometry, all other properties are essentially 3D p8 |. Due to the third condition 
the return probability is dominated by diffusive contributions [13|. When k-p = 0(a^ x ), 
Eq. (|D corresponds to 1 C !/d < Af < L/a, where M = kpS is the number of channels. 
Furthermore, in the metallic regime, one must have I C I C ^. Since £ ~ M£, this implies 

1<L/£<M. (9) 

In a finite system, the conditions (J8|) and (§) are not easily met simultaneously. We have 
performed exact diagonalizations for 128 x 4 x 4 and 128 x 8 x 8 lattices, using open boundary 
conditions (BC) in the longitudinal direction. In this case, P(x, x; 0) = 2X [1/3— x{L-x) j L?\. 

The results of our calculations are summarized in Figs. ^| and |3|. Figure |2| shows 
(5fp(E,x;t)) x in comparison with Eqs. and (|7|). We observe very good agreement. 

The value of P = (P(x,x;0)) x should be independent of (3. As can be seen in Fig. |2], the 
value of P does somewhat change with /3, albeit weakly. For narrower wires (128 x 4 x 4) 
we have observed that the ratio P\jPi (determined by fitting P = Pp independently for 



(3 = 1,2) becomes very small for small values of W (corresponding to X ~ 0.1) while it 
approaches unity for large values of W. A possible explanation for this deviation would 
be that for small M and small X, the condition (|8|) is no longer satisfied since t/a ~ M. 
Surprisingly, the form of the deviations is still very well described by Eq. ([/]) (not shown). 

Figure ^ shows the tails of fp(E,x;t) for weak disorder (X ~ 1) in comparison with 
Eqs. (|j), (||D and (|^). Since for very small values of X the tails decay so fast that we cannot 
reliably calculate them, we decreased the wire cross section and increased the value of W 
in Fig. |, thus increasing X. The quoted values of X were obtained by fitting Eqs. (||) and 
(||). The values thus determined differ somewhat between (3 = 1 and 2 (see Fig. ^) and this 
difference depends on the choice of E, rj and W. 

In summary we conclude that non- universal deviations from RMT statistics in quasi- ID 
wires are very well described by a NLSM not only in the body (see Fig. |2|) but notably also 
in the tails (see Fig. ||) of the distribution fp(E,x;t). 

2D case. In this case, according to Ref. |Tl| , corrections to /i are still given by Eq. (fTf), 
but now P = (P(r, r; 0)) r where P(r, r'\ iS) is the 2D diffusion propagator. For the tails of 



the distributions, the result of the NLSM is within a saddle-point approximation ||10| , |13 

/^r^expt-C^lnt) 2 ] (10) 

with 

C^(37rg/[A\n(L/£)]. (11) 

Note that according to ( |TTD the decay in the tails of Eq. ( |T0"D depends on (3, as in the quasi-lD 
case [Eq. (0)]. 

Recently, in Ref. Jl9| a different approach (a direct optimal fluctuation method [^(J) was 



used to calculate the tails of fp(E, r; t). According to Ref. [0, the tails of the distribution 
function are given by Eq. ([10]) but with Cp replaced by 

C = 7rg/[2\n(L/r )] (12) 
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[with ro = 0(kp 1 )] which differs in two respects from the prediction of the NLSM: First, £ 
in Cp is replaced by r in (|T2"D. Second, there is no /^-dependence. 

We have diagonalized the Hamiltonian ([!]) on a 100 x 100 lattice. Figure ^ shows correc- 
tions to fg for weak disorder. We find that the form of the deviations is very well described 
by Eq. (|7|). However, the values of Pp obtained for f3 — 1,2 differ by a factor < 1/2. A 



possible explanation for this deviation might be the following |Tj|: In the ballistic regime, P 
is no longer given by the diffusion propagator but may be dominated by a single-scattering 
expression which involves an additional factor (3/2 and thus Pi/P% < 1. It would be tempt- 
ing to deduce from this that in our case ballistic effects are important. However, this does 
not explain P\j Pi < 1/2. Numerical results |29| (albeit for rather small systems) indicate 
that in 3D, P\jP?, ~ 1/2 for the parameters chosen in |[29|| . 

Figure |5] shows the tails of the distribution functions. The tails are consistent with an 
exp[— C(lnt) 2 ] decay as predicted by Eq. (|T0|). We have thus verified that corrections to 
RMT distributions in 2D systems do give rise to log-normal tails. Our results suggest that 
the prefactor in the exponent does not depend on (3. This result is consistent with Eq. (|12"|). 
We have independently calculated the dimensionless conductance g using the usual linear 
response expression. We find that g is independent of f3 and g oc W~ 2 , as expected (inset of 
Fig. §). 

The inset of Fig. [| shows that C increases with decreasing disorder strength, as it should, 
albeit slower than W~ 2 . The increase of C for decreasing W is underestimated, because for 
weak disorder, the tails of the distributions have not reached the asymptotic regime. 

In summary we have reported on a study of rare events in disordered conductors, by diag- 
onalizing the tight-binding Hamiltonian (p]) and analyzing the probability of rare splashes of 
high wave-function amplitudes. Our ID results agree with those of 0. In the quasi-lD case, 
we have compared our data to an exact solution 0|| of the NLSM, and to a saddle-point 
approximation [1(J . We observe very good agreement between our results and those of the 
calculations based on the NLSM and thus conclude that the NLSM provides a quantitative 
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description of rare events in quasi- ID disordered conductors. In 2D systems, corrections to 
the body of the distribution functions are well described by results based on the NLSM, 
with a modified prefactor Pp. Moreover, we could verify that the tails of the distribution 
function in the vicinity of the metallic regime are log-normal. Thus our numerical inves- 
tigations, which are complementary to the analytical predictions, corroborate the overall 
picture suggested in Refs. [§-0 for quasi-lD and 2D systems. The coefficient describing 
the tails of wave-function distributions in 2D systems turns out to be independent of (3 for 
the parameters considered in this paper, as opposed to the quasi-lD case. This is consistent 
with the prediction of the direct optimal fluctuation method (Ref. ||19|| ). 
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FIGURES 
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FIG. 1. (f(E,x;t)) x for a chain of length L/a = 2000 with periodic boundary conditions, 
= and 7] = 0.2. The lines are determined from Eq. (||). The arrow indicates t c for L/£ = 4.76. 
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FIG. 2. (8fp(E,x;t)) x for a 128 x 8 x 8 lattice with E = -1.7, 77 = 0.01, and W = 1.0. For 
= 1 (o) and (3 = 2 (□), these results are compared to Eqs. (||),(||) ( ) and Eq. (|7|) ( ). 
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FIG. 3. f p {E, x\ t) for a 128 x 4 x 4 lattice; for X ~ 1, x ~ L/2, TV = 1.6, £ = -1.7, r/ = 0.01 
and /3 = 1, 2 compared to Eqs. @, (|) and (|). 
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Eq. (7), = 1, P = 0.058 



Eq. (7), = 2, P = 0.143 
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FIG. 4. (5fp(E, r; t)) r for a 100 x 100 lattice with periodic BC, for E = -1.6 and 77 = 0.005. 
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(In*) 2 

FIG. 5. Tails of (fp(E, r; t)) r for the same system as in Fig. f|, for W = 2, 3 and 4, for (3 = 1 
(o) and (3 = 2 (□). The lines show Eq. (11). The inset shows the fitted values of C versus W~ 2 , 
for p = 1 (A) and 2 (y) and the dimensionless conductance g for /3 = 1 (o) and 2 (□). The dashed 
line indicates W -2 behaviour. 
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